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ABSTRACT 



The southern magnetic Bp star HD 137509 exhibits complex rotational modulation of the longitudinal field and other magnetic observables. 
Interpretation of this magnetic variability in the framework of the low-order multipolar field models suggests a very strong quadrupolar compo- 
nent to dominate the surface field topology of HD 137509. 1 have examined the high-quality VLT/UVES spectra of HD 137509 and discovered 
resolved Zeeman split components in some of the spectral lines. The inferred mean surface field modulus, (5> = 29 kG, agrees with the multi- 
polar model predictions. This confirms the presence of an extremely strong non-dipolar magnetic field in HD 137509 and establishes this star 
as the object with the second-largest field among magnetic chemically peculiar stars. 
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1. Introduction 

Significant fraction of the middle and upper main sequence 
stars of the B-F spectral classes shows unusual surface proper- 
ties. The chemical composition of the outer layers in these stars 
departs significantly from solar, with elements distributed non- 
uniformly both horizontally over the stellar surface and verti- 
cally in the line forming region. Many of these chemically pe- 
culiar (CP) stars also possess a strong, well-ordered magnetic 
field. Unlike the rapidly evolving, complex spotted magnetic 
topologies of active late-type stars, CP-star magnetic fields are 
stable and cover the whole stellar surface. The field strength 
ranges from a few hundred G up to ^ 35 kG. The field is frozen 
in the stellar plasma and its observed changes are not intrinsic 
to the star but rather a result of the rotational modulation of the 
visible projection of the stellar magnetic field geometry. 

The structure and origin of the global magnetic fields of CP 
stars is a matter of intense debate. The early observational mag- 
netic field studies focused on the measurements and establish- 
ing the form of the rotational modulation of the disk-integrated 
line of sight field component - longitudinal field {B^). The vast 
majority of the magnetic CP stars exhibit smooth, single- wave 
variation of longitudinal field during the stellar rotation cycle. 
This was considered as an evidence for the dominant dipo- 
lar contribution to the field topology (e.g. Borra & Landstreet 
I1980> . However, after several additional magnetic observables 
were incorporated in the field modelling procedure (Landstreet 
& Mathys 2000 Bagnulo et al. 2002), it was found that purely 
dipolar fields are often insufficient to explain observations, 
hence more complex quadrupolar and octupolar fields had to 
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be introduced. Moreover, recent detailed analysis of the four 
Stokes parameter spectra of CP stars (Kochukhov et al. l2004> 
revealed that in some stars exceedingly complex magnetic ge- 
ometries may give rise to a deceptively simple, nearly sinu- 
soidal (Bz) variation. 

Thus, despite a certain success of the phenomenological 
low-order multipolar field models, it remains unclear whether 
the suggested quadrupolar and octupolar deviations from the 
basic dipolar structure reflect the physical reality or represent 
merely a mathematical approximation of the unresolved small- 
scale field. Solution of this problem is certainly of great in- 
terest because it will provide additional important constraints 
for the theoretical models of the field generation and evolu- 
tion in CP stars. I suggest that a better understanding of the 
nature of non-dipolar fields in early-type magnetic stars can be 
gained through detailed investigation of stars with the domi- 
nant non-dipolar field components. Up to now only two ob- 
jects with extraordinary double-wave longitudinal field curves, 
hence clearly non-dipolar fields, have been found. These stars 
are HD 37776 (Thompson & Landstreet 198J^, Khokhlova et 
al. l2000t and HD 137509 ( Mathy s & Lanz[l997). The third 
star, HD 133880 (Landstreet 19901 shows a strongly distorted 
single-wave {B^) variation and is often referred to as another 
example of the star with an unusually strong non-dipolar field. 

One could hope to constrain the structure of non-dipolar 
fields in these objects and distinguish a low-order multipolar 
field from a more complex, spot-like, magnetic topology with 
the help of the disk-integrated surface field strength {B) mea- 
surement. The low-order multipolar field scenario predicts that 
a relatively low-amplitude non-sinusoidal longitudinal field 
variability is produced by a quadrupolar or octupolar field of 
an enormous strength. For instance, a 4 kG peak-to-peak mod- 
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ulation of {B^) in HD 37776 (Thompson & Landstreet IT9831 
should be interpreted as a signature of > 50 kG surface field 
(Bohlender 1994 1. In contrast, discarding a low-order multipo- 
lar field parameterization, one can successfully reproduce the 
observed non-sinusoidal {B^) phase curves with a high-contrast 
surface field of a modest strength (only a factor of 1.5-2 larger 
than the {B^) extremum). 

In this paper I try to constrain different possible non-dipolar 
field structures in the Bp star HD 137509. Of the three early- 
type stars with distinctively non-dipolar magnetic field configu- 
rations, HD 137509 has the longest rotation period, lowest pro- 
jected rotational velocity and, therefore, presents the best op- 
portunity for the direct diagnostic of the surface field strength 
through the high-resolution spectroscopic observations and line 
profile analysis. 

The rest of the paper is organized as follows. In Sect. |2] 
I summarize the physical properties of HD 137509 and refine 
the value of the stellar rotation period. Sect. |3] overviews pre- 
vious magnetic field measurements and presents interpretation 
of the magnetic field variability in terms of the low-order mul- 
tipolar field topologies. Detection of the resolved Zeeman split 
lines in the spectrum of HD 137509 and the first field strength 
measurement are reported in Sect.|4] The study concludes with 
a summary and discussion in Sect.|5] 

2. Atmospheric parameters and rotation period 

HD 137509 (HIP 7601 1, NN Aps) was assig ned a B9p SiCrFe 
spectral classification by Cowley & Houk (I1975> . These au- 
thors mention strong lines of Si ii, Cr n, Fe ii, Ti ii, and the ab- 
sence of He I features. Atmospheric parameters of HD 137509, 
Teff = 1 2600 + 200 K and log g = 4.06 + 0. 1 1 , were determined 
by Hunger & Groote (1999'l using the IR flux method and com- 
parison with the stellar evolutionary tracks. On the other hand, 
Hubrig et al. ( l2000b suggested = 11900 + 600 K from 
the published spectral type and deduced log^ - 3.87 + 0.16 
from the Hipparcos parallax and theoretical stellar models. I 
find that the Wfi and Hy hydrogen lines in the UVES spectrum 
of HD 137509 discussed below (Sect. |4|i are best fitted with 
Jeff = 12750 + 500 K and \ogg = 3.8 + 0.1 if the enhanced 
metallicity [MjH] = -i-l.O atlas9 (Kurucz FT^^ stellar model 
atmosphere is employed. 

Photometric variability of HD 137509 was detected by 
Loden & Sundman J1989> . Mathys J1991> discovered a strong 
reversing longitudinal magnetic field and pointed to spectac- 
ular variations of the Si ii and iron-peak element line intensi- 
ties, confirming earlier tentative spectrum variability reports by 
Cowley & Houk (T975'i and Pedersen (7979*1. Mathys & Lanz 
J1997i) have combined all available magnetic measurements 
and obtained new time-series photometry in the Geneva sys- 
tem. Simultaneous consideration of these observables allowed 
the authors to obtain a refined value P - 4.4916 + 0.0002 d for 
the rotation period of HD 137509. 

In addition to the previous reports of periodic brightness 
modulation, HD 137509 also shows pronounced variability in 
the Hipparcos epoch photometry (ESA [T997b . The Hipparcos 
satellite observed the star over the period of 3.2 years, obtain- 
ing 192 individual photometric measurements. Fourier analysis 
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Fig. 1. Time-series analysis of the Hipparcos epoch photom- 
etry of HD 137509. Upper panel: amplitude spectrum of the 
photometric variation. The inset shows the window function 
of the Hipparcos observations. Lower panel: photometric mea- 
surements folded with the period P - 4.49257 d. The solid 
curve illustrates the least-squares fit by a sinusoid and its first 
harmonic. 

of these data is illustrated in Fig.^ The only statistically sig- 
nificant peak is detected at the frequency v = 0.223 d The 
least-squares fit by a cosine function and its first harmonic 

Vh^Vq + Vi cos [2n{T -Tq)/P + ^,] 
+ V2 cos [47t{T - Tq)/P + ip2] 

yields P = 4.49257 + 0.00029 d, Vq = 6.846 mag, V, = 0.019 ± 
0.001 mag, V2 - 0.012 ±0.001 mag. Here the epoch of positive 
magnetic extremum JD - 2448344.575 (Lanz & Mathvs lT99ll 
was adopted as the reference time Tq. 

The new estimate of the rotation period differs by < 3cr 
from the value derived by Mathys & Lanz J1997> and has 
similar formal precision. Moreover, the present analysis of the 
Hipparcos photometric time-series overcomes aliasing problem 
typical of the previous photometric observations (see fig. 1 in 
Mathys & Lanz 1 19971 and hence gives more confidence in the 
derived period. 

The brightness variation of HD 137509 as measured by 
Hipparcos correlates remarkably well with the magnetic field 
modulation. Comparison of the photometric curve (Fig.^ and 
(Bz) variability (Fig.|2j indicates that the star shows maximum 
brightness at the phase of positive magnetic extremum {tp - 
0.0). The secondary longitudinal field maximum at = 0.5 
also has its counterpart in the photometric curve. 
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3. Magnetic field topology 

The first longitudinal field measurements were obtained for 
HD 137509 by Mathys (1991 1 using metal line spectropo- 
larimetry. Additional (Z?z> observations based on the photopo- 
larimetric measurements in the line were acquired by 
Bohlender et al. ( 1993 1- Variation of other integral magnetic 
observables, the crossover Ve sin i{xBj) and the mean quadratic 
field {Bl + b2)1/2^ ^as studied by Mathys ( ll995al[T9^5bi l and 
Mathys &Hubrig( 19971. 

All longitudinal field measurements available for 
HD 137509 are illustrated in Fig. |2] In total there are 16 
(Bz) determinations, but only four of them correspond to a 
> 3cr field detection. Nevertheless, the longitudinal field vari- 
ation of HD 137509 is reasonably well defined and exhibits a 
remarkable double-wave shape, with a peak-to-peak amplitude 
of Si 3.5 kG and unusually narrow positive field extremum at 
rotation phases if - 0.8-1 .2. An even more pronounced double- 
wave variation is evident in the crossover measurements (see 
Fig-13- The mean quadratic field reaches up to 36 kG, which 
is the largest {B\ + B^y^^ found in magnetic CP stars (Mathys 
I1995bl l. These measurements provide a hint for the presence 
of a very strong surface field in HD 137509. However, the 
quadratic field diagnostic by itself is rather inconclusive 
because it relies on the interpretation of subtle magnetic 
line broadening effects and neglects horizontal chemical 
inhomogeneities. For several spotted magnetic CP stars very 
large quadratic fields inferred by Mathys ( 1995b I could not be 
confirmed by the subsequent detailed line profile analyses of 
high-resolution spectra (Kupka et al. 119961 Kochukhov et al. 
[20041 . 

Uncertainties of the magnetic observables notwithstanding, 
I have attempted to probe the field structure of HD 137509 in 
the framework of the low-order multipolar magnetic geome- 
tries and using all available measurements of (B^), Ve sin i{xB^), 
and (BI + B^y^^. The purpose of this exercise is not to study 
details of the stellar magnetic field topology, but to assess 
the range of the disk-averaged surface field strength {B} ex- 
pected for such simple models. The multipolar field modelling 
procedure applied to HD 137509 is similar to the analysis by 
Landstreet & Mathys 12000) and Bagnulo et al. J2002t . Having 
no a priori reason to prefer one or another form of the multi- 
polar field parameterization suggested in the literature, I have 
tested the three most widely used types of magnetic models. 
The first one (model A) is identical to the field parameteriza- 
tion used by Landstreet & Mathys (^000). The field is rep- 
resented by an aligned superposition of dipole, axisymmetric 
quadrupole and octupole. The 6 free model parameters include 
the respective field strengths B^, Bq, Bq, two angles, /3 and y, 
giving orientation of the magnetic field axis, and inclination / 
of the stellar rotation axis. The second model (B) disposes of 
the field symmetry assumption and approximates the magnetic 
structure by a sum of dipole and linear quadrupole, not nec- 
essarily aligned with respect to each other (Khokhlova et al. 
|2000l. This model has 7 free parameters: the field intensities 
Bd and Bq, four angles specifying orientation of the dipolar and 
quadrupolar components, and stellar inclination. Finally, model 
C is the one employed by Bagnulo et al. (i2002i . It is similar to 
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Table 1. Parameters of the best-fit multipolar magnetic field 
models of HD 137509. 



Parameter 


Model A 


Model B 


Model C 


i(°) 


81 + 15 


78 + 7 


64 + 21 


Bd (kG) 


3.1 ± 1.8 


4.9 ± 1.3 


5.5 + 2.2 


/?(") 


64± 12 


123 ± 26 


84+16 


rC) 


372 ±3 


367 ± 16 


372+ 11 


B, (kG) 


41.9 ± 1.3 


-46.8 + 1.3 


51.6 + 2.9 


/3i n 




43 ±6 


47 + 18 


Pin 




= A 


114+18 


n n 




81 ±4 


125 + 9 


72 D 




= 71 


216+ 10 


Bo (kG) 


-0.5 + 2.1 






xl 


1.36 


1.34 


1.16 
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Fig. 2. Observed variation of the magnetic field moments in 
HD 137509 and predictions of the best-fit multipolar mod- 
els of the magnetic field topology. The longitudinal field, 
quadratic field, and crossover observations are taken from 
Mathys ( .1991. 1995a L995bJ, Math ys & Hubrig ( .1997J {filled 
circles), and Bohlender et al. ( I1993I I (open triangles). A single 
data point shown in the lower right panel represents the field 
modulus measurement reported in the present paper. Smooth 
curves show the fit by multipolar models discussed in the text: 
model A {solid line), model B {dashed line), and model C 
{dash-dotted line). 

model B, but includes non-linear quadrupole, whose geometri- 
cal properties are now defined with the 4 angles. Since the field 
in HD 137509 is dominated by a higher-order magnetic com- 
ponent, it is convenient to reckon the azimuth angles yi and 72 
from the same zero-phase plane as the dipolar azimuth angle y. 
The projected rotational velocity, Ve sin /, influences the ampli- 
tude of the Ve sin i{xB^) variation and is adjusted for each of the 
three magnetic models to fit the crossover observations. 

Parameters of all magnetic models were optimized using 
the least-squares procedure. Resulting best-fit multipolar field 
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Fig. 3. Comparison of the observed (thin line) and computed (thick line) spectrum of HD 1 37509. Synthetic spectrum is computed 
for VeSin/ = 20 kms ' and {B}-29 kG using the synthmag code. Zeeman patterns are schematically illustrated for the Sin 
6347.1 1, 6371.37 A and Fe ii 6357.16 A lines. The length of vertical bars is proportional to the relative strength of the respective 
Zeeman components; n components are shown above the horizontal dotted line, and cr-t components below this line. 



characteristics are listed in Table ^ Theoretical phase curves 
of the four magnetic observables are compared with available 
observations in Fig. |2 The three magnetic models yield re- 
duced of 1.16-1.36, indicating acceptable fitting of the data. 
The general superposition of dipole and non-linear quadrupole 
(model C) provides a marginally better description of the mag- 
netic measurements. It is encouraging to find consistency in 
some of the geometrical parameters recovered in the fit with 
different multipolar parameterizations. All models favour large 
inclination angle, / = 64°-81°, and suggest similar strength and 
orientation for the dipolar component. The strength of the dom- 
inant quadrupolar component ranges from 42 to 52 kG - at least 
9 times larger than the dipolar contribution. The projected ro- 
tational velocity deduced from the crossover observations is in 
the interval 9-13 kms which is smaller than the spectroscop- 
ically determined Vg sin / (see Sect. 1^. The stellar radius esti- 
mated using the aforementioned range of possible stellar incli- 
nation and the rotation period of HD 137509 is also unrealisti- 
cally small (R x R^) for Vg sin / » 10 km s ' . On the other hand, 
one finds R » 1.8-3.0/?© if Ve sin / = 20-30 km s ' is adopted. 

The three magnetic models predict nearly identical rota- 
tional modulation of (Bz) and Vg sin i{xB^}. Expected variation 
of (BI + B^y^^ is marginally different, but corresponding the- 
oretical curves cannot be distinguished based on observations 
because of the large errors of the quadratic field measurements. 
Thus, no unique multipolar field model can be derived for 
HD 137509 with the available magnetic observations. This situ- 
ation is often encountered in the analysis of the disk-integrated 
magnetic observables of magnetic CP stars when using some 
specific parameterization of magnetic geometry (Bagnulo et al. 
12002 '^ or when comparing results obtained with different multi- 
polar parameterizations (Kochukhov .2004j . At the same time. 



one can see that the mean surface field variation expected for 
the alternative magnetic models of HD 137509 is substantially 
different, especially for the rotation phase interval (p =0.6-1.0. 
Predicted (B) ranges from 22 to 3 1 kG. CoiTesponding Zeeman 
splitting of the magnetically sensitive spectral lines should be 
comparable or larger than the rotational Doppler broadening 
of the HD 137509 spectrum (vg sin / = 28 + 5km s"' according 
to the estimate by Mathys ( 1995b|). Thus, if one believes that 
the low-order multipolar field models provide an adequate de- 
scription of the non-dipolar field in HD 137509, there are good 
prospects for detecting Zeeman split lines in high-resolution 
optical spectra of this star. 

4. Zeeman split lines in the spectrum of HD 137509 

High-resolution spectra of HD 137509 were obtained with the 
UVES instrument of the ESQ VLT on July 12, 2001 (JD = 
2452102.60, phase 0.67 for the 4.4916 d rotation period). 
Observations were carried out in the context of the UVES 
Paranal Observatory Project^ (ESO DDT program 266. D- 
5655, Bagnulo et al. I2003t . The spectra were recorded in 
dichroic modes, using the slit width of 0.5". A spectral resolu- 
tion of about 80000 and almost complete spectral coverage of 
the wavelength interval from 3070 to 10400 A were achieved. 
The signal-to-noise ratio of the HD 137509 spectrum is above 
300 for the optical region. Reduction of the UVESpop obser- 
vations is described by Bagnulo et al. (2003). 

Inspection of the UVES spectra of HD 137509 revealed a 
significant distortion in the profiles of many spectral lines. In 
addition to the overall asymmetry of many features, some lines 

' jhttp : //www ■ eso . org/uvespop/| 



O. Kochukhov: Remarkable non-dipolar magnetic field of the Bp star HD 137509 



5 



exhibit resolved absorption components which are consider- 
ably narrower than expected for the normal late B stellar spec- 
trum with a v'e sin / «30 km s ' Doppler broadening. This may 
be a signature of the Zeeman splitting in a very strong magnetic 
field or a manifestation of an inhomogeneous surface chemical 
distribution. To separate these effects, I have investigated sev- 
eral groups of spectral lines, belonging to the same ions, but 
with different sensitivity to magnetic field. 

Remarkable proof of the presence of a very strong magnetic 
field in HD 137509 comes from the analysis of the red Si ii dou- 
blet at A 6347 and 6371 A. These lines originate from the same 
lower atomic level and have similar strength. Nevertheless, 
as evident from Fig. |3l profile shape of the two Si ii lines is 
very different. The 6371 A line shows an asymmetric double 
structure, whereas the Sin 6347 A lines exhibits no apparent 
anomaly. Obviously, horizontal abundance spots cannot be re- 
sponsible for such a drastically different shape of the two simi- 
lar lines. Consideration of the Zeeman effect allows one to un- 
derstand the puzzling behaviour of the Si ii doublet. The Si ii 
6371 A line has a relatively simple, nearly doublet, Zeeman 
splitting pattern with two n and two cr components. On the 
other hand, the 6347 A line splits into 6 close Zeeman com- 
ponents, which are blended together and do not produce re- 
solved structure in HD 137509. In Fig.|3ll compare observa- 
tions and theoretical spectrum calculated with the magnetic 
spectrum synthesis code synthmag (Piskunov I1999t . I have 
adopted T^ff = 12750 K, log^ - 3.8, Vg sin/ = 20 kms"', and 
(B) - 29 kG. The Zeeman splitting of Si ii 637 1 A line and 
broadening of the Si n 6347 A and Fe ii 6357 A lines are re- 
produced fairly well, although a homogeneous magnetic field 
adopted in the synthmag calculations is insufficient to model 
relative strengths of the Sin 6371 A Zeeman components. 

A number of doubly ionized Nd lines can be identified in 
the spectrum of HD 137509. Some of these features show con- 
spicuous narrow components in the triplet-like splitting pat- 
terns (Fig.E). Detailed investigation of several lines confirms 
this to be a magnetic field effect. The Ndni lines with sig- 
nificant Lande factors and quasi-triplet Zeeman splitting (e.g. 
5845 and 6145 A) exhibit resolved components. A well-defined 
triplet structure is also observed for an unidentified line at 
6273.67 A. According to the Nd laboratory line list compiled 
by H. Crosswhite, this line belongs to the Nd iii spectrum (T. 
Ryabchikova, private communication). On the other hand, as 
expected from its small Lande factor and complex Zeeman pat- 
tern, the resonance Nd iii 6550 A line shows no magnetic split- 
ting. Resolved Zeeman components of the Ndiii lines suggest 
Ve sin / ^ 10 km s ' , which is substantially smaller than the pro- 
jected rotational velocity inferred from the Si ii and other spec- 
tral lines. This anomaly may be attributed to a high-contrast 
inhomogeneous horizontal distribution of Nd. In such a situa- 
tion Nd lines will sample a fraction of the visible stellar surface 
and, therefore, will be less broadened by rotation. 

To improve the accuracy of the magnetic field strength mea- 
surement, I have determined (B) from several lines showing 
Zeeman components: Sin 6371, Feii 5018, Ban 6141, Ndm 
5845 and 6145 A. Fitting superposition of Gaussian profiles to 
these five lines yields (B) = 28.7 + 0.8 kG (the quoted uncer- 
tainty is the standard error of the mean for the measurements 
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Fig. 4. Zeeman resolved Ndm lines in the spectrum of 
HD 137509. The four panels show lines with different Zeeman 
splitting patterns. The strength and position of the n and cr 
Zeeman components for {B) - 29 kG are illustrated as in 
Fig E] The classical triplet pattern with z - 0.9 is adopted for 
the unclassified line Ndm 6273.67 A. 

obtained from individual lines). The aforementioned spectrum 
synthesis modelling of the Si n doublet is in agreement with 
this estimate of the surface field strength. 

The first field modulus measurement obtained for 
HD 137509 in the present paper is compared with the multi- 
polar field model predictions in Fig.|2l The maximum attained 
field strength is consistent with observations for all three mod- 
els. The non-axisymmetric superposition of the dipole and lin- 
ear quadrupole (model B, see Sect.|3} shows the best agreement 
with the observational estimate of (B), whereas the axisymmet- 
ric model A is probably ruled out. 

Magnetic spectrum synthesis with the synthmag code was 
used to assess chemical abundance pattern of HD 137509. Even 
when the very strong magnetic field is taken into account, one 
finds a large enhancement of the Si and iron-peak element 
concentration with respect to the solar chemical composition: 
[Si/H] = -hO.8, [Fe/H] = -h1.4, [Cr/H] = -h1.9, and [Ti/H] = -h2.6. 
In contrast, several light elements are strongly underabundant: 
[He/H] < -3.5, [Ca/H] = -2.2, [Mg/H] = -1.2. 

5. Discussion 

In this paper I have presented analysis of the magnetic field 
in the unique peculiar B-type star HD 137509. This object is 
one of the very few early-type stars with distinctively non- 
dipolar magnetic field geometry. Interpretation of the previ- 
ous magnetic measurements of HD 137509 with the help of 
the low-order multipolar field models suggests that the dom- 
inant contribution to the surface field in this star comes from 
the > 40 kG quadrupolar component. Mutipolar models predict 
the disk-averaged surface field strength in the range from 22 
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to 3 1 kG. My investigation of the high-resolution UVES spec- 
tra of HD 1 37509 revealed resolved Zeeman splitting in several 
spectral lines. The inferred field strength of 28.7 kG is in accord 
with the multipolar model predictions. 

Unlike several other cases when a large quadratic field de- 
duced for spotted magnetic CP stars contradicts high-resolution 
line profile studies (see Kupka et al. (1996) for a Cir and 
Kochukhov et al. ( 2004 ) for HR 3831), the quadratic field mea- 
surements of HD 137509 appear to be more reliable. This is 
likely to be a result of a comparatively strong influence of the 
magnetic field on the line profile shapes in HD 137509. In stars 
with weaker fields effects of the horizontal chemical inhomo- 
geneities often dominate the line profiles, which distorts the 
quadratic field diagnostic (Kochukhov 120041 ) and sometimes 
leads to a spurious detection of strong magnetic fields. 

The present study establishes HD 137509 as an outstand- 
ing object with the second-largest field among CP and main 
sequence stars in general. The only star with a stronger field is 
the well-known Babcock's star (HD 215441). The latter shows 
varying field modulus between 32 and 34 kG (Preston 1969). 
However, orientation of the rotation axis and magnetic field 
geometry of HD 215441 are not favourable for detailed anal- 
ysis of the field topology (Landstreet et al. 11989 ). In con- 
trast, large inclination angle and unusual, very strong field of 
HD 137509 open remarkable possibilities for the detailed mag- 
netic field mapping (Piskunov & Kochukhov l2002t and for the 
analysis of relation between chemical spots and magnetic struc- 
tures. Outstanding chemical enrichment of the HD 137509 at- 
mosphere allows to use this star as a test ground for the appli- 
cation of the new generation model atmospheres of early-type 
stars, incorporating effects of anomalous chemical composition 
and strong magnetic field (Shulyak et al. .2004. Valyavin et al. 
l2004l Kochukhov et al. 2005 ). 

The measurement of the field strength in HD 137509 is the 
first case when we are able to impose a strong constraint on 
the fundamental field characteristic of an early-type star with 
a non-dipolar field topology. Good agreement between (B) de- 
termination and the outcome of the multipolar field analysis 
supports the common (through not as yet thoroughly verified) 
assumption that deviations from the basic dipolar field structure 
in the early-type magnetic stars has the form of quadrupole or 
octupole. Justification of the mutipolar field modelling frame- 
work hints that another star with the complex, double-wave 
variability of the longitudinal field - HD 37776 (Thompson 
& Landstreet 1985) - also hosts a similarly large magnetic 
field, or may even possess the strongest field of all known non- 
degenerate stars. Unfortunately, the rapid rotation of HD 37776 
precludes us from detection of the Zeeman splitting in spectral 
lines and hence leaves no chance of direct confirmation of the 
presence of an extraordinarily large field in this star. 
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